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Scattering probe particles from a quantum system can provide experimental access to informa- 
tion about the system's state. However, measurement backaction and momentum transfer during 
scattering changes the state of the system, potentially destroying the state of the system we wish to 
probe. Here we investigate how to probe the system's initial state even in the presence of backaction 
and momentum transfer. We show that summing the scattering distributions of an ensemble of mea- 
surements reveals the initial state scattering pattern even when each of the individual measurements 
completely destroys the initial state. This procedure is effective provided the scattering takes place 
on a timescale that is short compared with the free evolution of the system. 

PACS numbers: 03.65.Ta 37.10.Jk 03.65.Nk 



I. INTRODUCTION 

Scattering of particles from matter is a fundamental 
method used to reveal a system's properties. Ruther- 
ford's revolutionary discovery of the atomic nucleus in- 
volved scattering alpha particles from gold foil [1], while 
scattering of x-ray photons has led to numerous discov- 
eries about crystalline structure, such as the structure of 
DNA [2]. The range of phenomena investigated by scat- 
tering is broad, where for example neutron scattering has 
been used to probe the nature of exotic forms of matter 
such as superfluid helium and high temperature super- 
conductors [3] , and scattering of optical photons has been 
used to probe correlations within Bose-Einstein conden- 
sates [4]. In all these cases probe particles scatter from 
the matter because of density or charge fluctuations in 
the system or by creating or destroying excitations in 
the system. Measuring the scattering pattern then yields 
valuable information about the system's spatial structure 
and excitations, and furthermore interference between 
different scattering channels reveals information about 
correlations within the system. 

Typically in experiments and theoretical discussions 
involving scattering, the effect of the probe on the sys- 
tem is assumed to be weak and the scattering pattern 
is calculated in the Born approximation from the ini- 
tial system state [5]. However during each scattering 
event energy and momentum may be transferred from 
the probe to the system. Furthermore the act of measur- 
ing the probe's scattering pattern may produce quantum 
backaction on the state of the system [6, 7]. These effects 
can combine to significantly alter the state of the system, 
indeed for two particles, scattering of light leads to local- 
ization of relative position while the particles' momen- 
tum states become entangled [8-10]. This complicates 
the use of scattering to probe the state of a system, as 
after a small number of scattering events the system may 
be in a completely different state from which it started. 

Here we model the evolution of manybody states 



caused by consecutive scattering and measurement, and 
see how quantum backaction changes the scattering dis- 
tribution for multiple particle systems. As we have previ- 
ously shown, scattering leads to localization of the system 
particles in position space and, in some regimes, super- 
positions of localized particle configurations [11]. Here 
we show that the average scattering pattern for an en- 
semble of experiments gives the pattern predicted for the 
initial state even if this initial state is completely de- 
stroyed in each individual scattering experiment. This 
is possible when the scattering measurements all occur 
on a timescale short compared with the free evolution of 
the system. Experimental access to the initial state of 
the system is hence possible even for systems that are 
easily perturbed. This is an important result that we ex- 
pect can be utilized, for example, in the field of quantum 
simulation, where experiments using cold atoms are used 
to simulate many-body physics [12, 13]. The outcome 
of each cold atom experiment can be measured by scat- 
tering light or matter from the system [14-20] and then 
the state of the system can be reproduced by combining 
results from multiple runs of the experiment. 



II. CONSECUTIVE SCATTERING FROM 
MANYBODY STATES 

The model of consecutive scattering that we use was 
presented in Ref. [11] and we reproduce the relevant re- 
sults here for completeness. We consider coherent scat- 
tering of a probe beam from a system of identical par- 
ticles, where the probe particles are assumed to be in 
approximately plain wave form with initial wave-vector 
ki. Through scattering the probe can acquire a new 
wave- vector k^, which is then associated with a mo- 
mentum transfer to the system of /ik = /i(ki — k^). 
For a weak probe-system interaction the probability of 
multiple scattering is low and we can then assume this 
momentum transfer occurs in a single scattering event. 
For an iV-particlc state ~ J rfRV'(R)|R), where 
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R = {ri,...,r7v} gives the coordinates of the N par- 
ticles, the scattering interaction resuhs in a combined 
probe-system state of 



fiRdk/S'(k,R)V'(R)|R) ® |k/) (1) 



where 



S'(k,R) 



for 
for 



k/ ^k. 



(2) 



This state is superposition of ah possible single scatter- 
ing events along with the possibility that the probe did 
not scatter and remains in |ki) [11]. The factor g is given 
by the strength of the probe-system interaction and for 
simplicity we assume g is independent of ky , which is the 
case when energy conservation restricts scattering so that 
|k/| ^ |ki| [5]. This occurs when the probe particles are 
much lighter than the system particles, as in the case of 
light scattering, or can occur as we see below in the case 
of cold atoms in optical lattices where the system par- 
ticles are confined to energy bands. We further assume 
the probe-system interaction is isotropic, as is the case 
for s-wave scattering in cold atoms. The treatment can 
be generalized to take into account anisotropic scatter- 
ing, such as the dipole pattern in light scattering, but this 
does not affect the qualitative behavior discussed below. 

Measurement of the final probe wave- vector, for ex- 
ample by imaging in the far field, causes a backac- 
tion on the system state collapsing this superposition 
IV') — > J dRS'(k,R)V'(R)|R). If a probe particle is 
observed to be scattered with wave- vector k/ then the 
new system state is a superposition of each system parti- 
cle having received momentum Kk from the probe, given 
the new state a broader relative momentum distribution. 
The superposition occurs as detecting the probe's final 
wave-vector reveals no information about which particle 
scattered it. As the system's state broadens in relative 
momentum space it localizes in relative position space 
[8]. However, because the scattering probability distri- 
bution is completely determined by the relative positions 
of pairs of system particles, the backaction of the mea- 
surement on the state will then preserve superpositions 
of states that have the same set of relative position vec- 
tors. Measuring the scattering distribution can then lead 
to spatial superpositions of the system particles as we 
describe in Ref. [11]. 

Due to the small value of g it is often the case that the 
probe particles do no scatter from the system. Naively 
one would assume that this does not affect the system, 
but the measurement of a non-scattering event also re- 
veals information about the state and leads to measure- 
ment backaction [8]. This is because some system con- 
figurations scatter less than others and observing non- 
scattering projects the state toward these configurations. 
The non-scattering amplitude ^(R) is fixed by the con- 
dition that probe particles must either be scattered or 
non-scattered, that is 1 = ^(R)^ -I- dk|S'(k, R)p = 
/dk|5(k,R)|2. 



Together consecutive scattering and non-scattering 
events lead to a dynamic evolution of the many-body 
state. Neglecting, for the moment, the free evolution of 
the system between each scattering event, the net result 
of the scattering evolution is that each scattering mea- 
surement provides a partial measurement of the relative 
positions of the system particles. In the limit of an in- 
finite number of scattering events the relative positions 
of the system particles are determined exactly, at least 
in principle, and the resulting state is an eigenstate of 
relative position, such as the completely localized state 
|R). The probability of the state ending in |R) is given 
by the initial state wavefunction |?/;(R)P and the scat- 
tering pattern then observed is |5(k, R)p. If we then 
do multiple repetitions of the experiment and image the 
scattering pattern for each, we will get an ensemble of 
scattering patterings corresponding to various |R) with 
probability |?/'(R)p. Summing the scattering patterns 
over the ensemble then gives the same result as the scat- 
tering pattern for the initial state (without backaction) 

/ dR|S•(k,R)?/'(R)|^ 

More generally for a finite number of 
scattering events, m, the state becomes 
/dR5(ki,R)...5'(k™,R)V'(R)|R). Then in 

an ensemble of measurements with m events 
each, any scattering event we choose, say 
J, has a probability distribution given by 
/ dki . . . dk,_idk,+i . . . dk,„|5(ki, R) . . . 5(k,„, R)VXR)P 
/ (iR|S'(kj, R)?/'(R)p, that is the initial state pattern. 
Combining the scattering patterns from an ensemble of 
measurements then gives an estimate of the initial state 
scattering pattern with a sampling error that scales with 
1/v^ where n is the number of experimental runs. This 
occurs despite the fact that scattering leads the initial 
state to be destroyed in each individual run. 



III. SCATTERING AND FREE EVOLUTION IN 
THE BOSE-HUBBARD MODEL 

The analysis above assumes the system does not evolve 
while the scattering measurements are taking place. To 
examine how the system's free evolution affects the mea- 
surement we now examine a specific example that dis- 
plays the key characteristics of the dynamic scatter- 
ing process. We consider scattering from a system of 
Bosons on a one-dimensional lattice described by the 
Bose-Hubbard model. This can be achieved in experi- 
ment by matter-wave scattering from a one-dimensional 
lattice of bosonic atoms in an optical lattice [16, 17, 21]. 
We take the lattice to have M sites and to be oriented 
along the y-axis, while the initial wave- vector of the probe 
is in the x-direction. For simplicity we only consider scat- 
tering within the a;y-plane, where the full three dimen- 
sional scattering is a straightforward generalization. 

At low temperature the atoms all reside in the lowest 
band of the lattice and the energy of the probe can be 
arranged so that excitation of the atoms in the lattice to 
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higher bands is neghgible [16]. The state of the system 
can then be described in terms of lowest band Wannier 
functions w{r—rj), where rj is the position of the jth lat- 
tice site [22]. When this condition is not met we expect 
to get similar results by taking into account the contri- 
butions from higher bands [18]. We assume that the lat- 
tice potential is strong enough that the overlap between 
neighboring Wannier functions is negligible. The state of 
the system can be expressed in terms of a number basis 
|n„) = \{nj,j — 1, . . . , M}), where is the number of 
atoms at site j and u uniquely identifies each basis state. 
After the mth scattering event we expand the state as 
I*™) = ■^/'tTl'^M)- Scattering from this state then oc- 
curs at angle 9 to the x-axis with probability 



(3) 



where 1(6) = / dre*''(«) "-|w(r)|2 and k(6') = fco(l - 
cos 9, — sin 9). 

Following a detection at 9 the many-body state is pro- 
jected into the new state 



1^ 



m+l/ 



1 

T 



MS), 



>u), (4) 



where Af normalizes the state. We note that the number 
basis states are eigenstates of this projection, and the 
scattering process will preserve any state that begins in 
a basis state. Moreover, some of the basis states produce 
the same scattering pattern as the relative positions of 
pairs of atoms in the lattice are the same, for example 
in the iV = M = 3 case, |201) and |102) resuU in the 
same light scattering. Superposition of these states are 
partially preserved by the projection, in that the weight 
of each state in the superposition remains the same after 
scattering but the phase relationship is changed. 

Alternatively detection events where the probe is not 
scattered occur with probability 



E 



(5) 



and 



■Au — 



\ 



1-1-/ M 



(6) 



Detecting a non-scattering event projects the state into 
the new state 



(7) 



which favors states with a lower probability of scattering. 

The dynamic scattering process can now be simulated 
using a quantum jump procedure [23, 24]. Taking the ini- 
tial state, we calculate the probability distributions for 



scattering and non-scattering. A pseudo-random num- 
ber is then used to determine if scattering occurs and in 
which direction. If it does then the projection in Eq. (4) 
is applied, if it does not the non-scattering projection in 
Eq. (7) is applied. In either case the many-body state is 
then normalized and evolves according to Bose-Hubbard 
Hamiltonian before becoming the input state and the 
process repeats. For simplicity we assume the scatter- 
ing events are evenly spaced by time dt. 

As an example of the dynamic process we look at the 
simple case of a three site lattice containing three atoms. 
In this case there are only ten basis states making it 
straightforward to track the development of the many- 
body state. In Figure 1 we show realizations of the dy- 
namic scattering process for the three site lattice in the 
superfluid regime of the Bose-Hubbard model. For this 
example we have set the coupling constant g = 0.1 and 
as a result the vast majority of detection events result 
from non-scattering, shown as = in Figure 1(a). In 
Figure 1(c) we show the development of the many-body 
state for the case of no system evolution between each 
scattering event, dt = 0. In this example the many-body 
state progresses toward a superposition of the states |201) 
and 1 102), two states which produce the same scattering 
pattern. As discussed above, continued scattering from 
this end superposition does not change the constituent 
basis vectors but does change the phase of the superpo- 
sition. In Figure 1(b) we see that detection at non-zero 
angle quickly reduces the overlap of the many-body state 
with the original ground state. We see that the overlap 
makes quantum jumps when a scattering event occurs 
and gradually evolves due to non-scattering events. 

Figure 1(d) and (e) demonstrate the effect of system 
evolution between each scattering event. For the exam- 
ple evolution in Figure 1(d) with dt = O.OOlJ/h the sys- 
tem progresses toward a superposition like that in Fig- 
ure 1(b), however after some time the system evolution 
causes the system to jump to another set of states. In 
Figure 1(e) with dt ~ O.OIJ/H the evolution is completely 
different with the system exploring a range of basis states 
over the simulation. This is detrimental to determining 
the initial state scattering pattern and, as we will see be- 
low, to do so we must limit the time over which scattering 
occurs to less than the natural system evolution time. 



IV. SCATTERING DISTRIBUTIONS FROM 
ENSEMBLE MEASUREMENTS 

By repeated simulation of the dynamic scattering pro- 
cess with dt = 0, we find that the simulations all set- 
tle into a final state after a small number of scatter- 
ing events, as in Figure 1(c). The end states are al- 
ways superpositions of eigenstates of the scattering pro- 
jection that produce the same scattering pattern. We 
also see numerically that the proportion of times a sim- 
ulation ends in a particular state is determined by the 
initial state. For example for initial state |4'o) we get the 
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FIG. 2. Comparison of scattering distribution summed over 
multiple simulations for a lattice of 9 sites containing 9 atoms 
for two different ratios of the Bose-Hubbard parameters U 
and J. (points) Number of event detections Nd{0) after (a) 10 
simulations of 10000 detection events each and (b) 1000 sim- 
ulations with 1000 detection events each, where the angular 
range is divided into 600 bins for event counting, (lines) Scat- 
tering distributions predicted from the initial ground states. 
Parameters used are g — 0.1 and fco — Tr/a. 
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FIG. 1. Development of an atomic state caused by 1000 scat- 
tering events for a lattice with three sites and three atoms, (a) 
Detected events, (b) Overlap of the many-body state with the 
initial ground state |^'o). (c)-(e) Modulus squared of the basis 
coefficients of the state |*™) for (c) dt = (b) dt = 0.001 J/ft 
and (c) dt — O.OlJ/fi. Parameters used are U/J = 0.05, 
g = 0.1 and fco = n/a, where a is the lattice site separation 
and U and J are the parameters of the Bose-Hubbard model. 



end state in Figure 1(c) with probability |(201|*o>P + 
|(102|«'o)p. 

Because the final state proportions are the same as the 
initial basis state probabilities, we can obtain the initial 
state scattering distribution if we average over an ensem- 
ble of scattering experiments. This occurs even though 
the initial state is completely changed in the scattering 
process. For m scattering events per run of the experi- 
ment and n runs, each individual run will push the state 
towards a particular set of basis vectors with a probabil- 
ity determined by the original state. Repetitions of the 
experiment then lead to a statistical determination of the 
initial state scattering distribution as we probabilistically 
sample the scattering distributions for the basis states 
that make up the initial state. To get the best estimate of 
the initial state pattern with limited recourses, if we are 
restricted only by the total number of scattering events 
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mn then the optimal method would be to take m = 1 in 
which case back-action does not play a part. In many in- 
teresting experiments however, doing an experiment run 
is more restrictive than scattering. Then the error in es- 
timating the distribution by sampling n runs scales with 
1 / sjn. Furthermore, when ri 3> 1 increasing the number 
of scattering events per run increases accuracy with the 
same scaling in m and the accuracy of determining the 
photon distribution scales with 1 / y/mn. 

In Figure 2 we show how summing over multiple sim- 
ulations reproduces the initial ground state pattern for 
a lattice with 9 sites and 9 atoms. Despite there being 
6420 unique scattering patterns from the associated basis 
vectors, sampling only 10 is sufficient to distinguish the 
two different ground states of the Bose-Hubbard model 
we have show here. In Figure 2(a) each of the 10 sim- 
ulations includes 10000 events of which 10% (5%) were 
scattering events for U / J = 0.1(10). In Figure 2(b) we 
show the distribution from summing 1000 simulations 
with 1000 events each, in which case the error in the 
scattering pattern compared with the exact initial state 
pattern is reduced by a factor of approximately vTO in 
accordance with the 1/ ^mn scaling. 

The above results applied with no system evolution be- 
tween scattering events. As we saw in Figure 1(e) system 
evolution causes the state to eventually progress through 
the basis states rather than just settling into one group. 
One might wonder whether the basis states might then 
be sampled in proportion to the initial state and allow us 
to get the ground state scattering pattern without hav- 
ing to re-run the experiment multiple times. This is not 
the case and in fact for simulations long compared to 
the system's natural evolution time the basis states end 
up being sampled equally and the scattering pattern is 
independent of the initial state as shown in Figure 3(a). 
There are intermediate regimes however where scattering 
still typically results in a single group of basis states and 
the ground state can be recovered by summing over ex- 
periments. Figure 3(b) shows the result of summing over 
1000 simulations where each simulation lasted at total 
time of one natural unit of system time. In this case the 
summed scattering distribution closely matches the ini- 
tial state pattern. For the same parameters but a total 
time of 10 units, Figure 3(c), we see a significant depar- 
ture from the initial state pattern and the patterns from 
the two states begin to merge. 

V. CONCLUSION 

In conclusion, we have shown how scattering in the 
presence of measurement backaction scattering leads the 
initial state of the system to be completely destroyed. 
When the system does not have time to evolve this leads 



the system towards eigcnstates of the scattering measure- 
ment. Each individual experiment then yields a scatter- 
ing pattern corresponding to a particular eigenstate with 
probabilities determined by the initial state. Summing 




FIG. 3. Comparison of scattering distributions with different 
time dt between scattering events for a lattice of 3 sites con- 
taining 3 atoms, (points) Number of event detections Nd{0) 
for (a) one simulation of 1000000 events with dt = IJ/h, 
(b-c) 1000 simulations of 1000 detection events each for (b) 
dt — O.OOlJ/ft, (c) dt — O.OlJ/ft, where the angular range 
is divided into 600 bins for event counting, (lines) Scattering 
distributions predicted from the initial ground states. Param- 
eters used are g = 0.3 and ko = n/a. 



these scattering patterns then gives an estimate of the 
initial state scattering pattern. This allows the initial 
state to be probed even though it is completely destroyed 
in each run of the experiment. We have further shown 
that when the scattering happens over a time that is 
larger than the natural evolution of system this method 
no longer yields the initial state pattern, and thus to op- 
timize measurement of the initial state we must scatter 
all particles within the natural evolution time or less. 
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